On Optimal Short Recurrences for Generating Orthogonal Krylov Subspace Bases
نویسندگان
چکیده
In this talk I will discuss necessary and sufficient conditions on a nonsingular matrix A, such that for any initial vector r0, an orthogonal basis of the Krylov subspaces Kn(A, r0) is generated by a short recurrence. Orthogonality here is meant with respect to some unspecified positive definite inner product. This question is closely related to the question of existence of optimal Krylov subspace solvers for linear algebraic systems, where optimal means the smallest possible error in the norm induced by the given inner product. The conditions on A were first derived and characterized more than 20 years ago by Vance Faber and Tom Manteuffel (SIAM J. Numer. Anal., 21 (1984), pp. 352– 362). Their main theorem is often quoted and appears to be widely known. Its details and underlying concepts, however, are quite intricate, with some subtleties not covered in the literature. The talk will be based on joined work with Zdeněk Strakoš [1], and with Vance Faber and Petr Tichý [2]. Acknowledgement: The work was supported by the Emmy Noether Pro-gram of the Deutsche Forschungsgemeinschaft (J. Liesen and P. Tichý), and bythe Czech National Program of Research “Information Society” under project1ET400300415 (Z. Strakoš and P. Tichý) and the Institutional Research PlanAV0Z10300504 (Z. Strakoš). References[1] J. Liesen and Z. Strakoš, On optimal short recurrences for generating orthogonalKrylov subspace bases, to appear in SIAM Rev.[2] V. Faber, J. Liesen and P. Tichý. The Faber-Manteuffel Theorem for linear operators,submitted to SIAM J. Numer. Anal.
منابع مشابه
The Faber-Manteuffel Theorem for Linear Operators
A short recurrence for orthogonalizing Krylov subspace bases for a matrix A exists if and only if the adjoint of A is a low degree polynomial in A (i.e. A is normal of low degree). In the area of iterative methods, this result is known as the Faber-Manteuffel Theorem (V. Faber and T. Manteuffel, SIAM J. Numer. Anal., 21 (1984), pp. 352–362). Motivated by the description in (J. Liesen and Z. Str...
متن کاملOn Short Recurrences in Optimal Krylov Subspace Solvers:
To solve large sparse linear systems of equations computational fast methods are preferable above methods such as GMRES. By exploiting structure of the involved matrix one can in certain cases use short recurrences to create an optimal Krylov subspace method. In this report we explain some of the theory regarding short-recurrence methods for solving linear systems. A central object in this theo...
متن کاملMultiple Recurrences and the Associated Matrix Structures Stemming From Normal Matrices
There are many classical results in which orthogonal vectors stemming from Krylov subspaces are linked to short recurrence relations, e.g., three-terms recurrences for Hermitian and short rational recurrences for unitary matrices. These recurrence coefficients can be captured in a Hessenberg matrix, whose structure reflects the relation between the spectrum of the original matrix and the recurr...
متن کاملComputing Approximate Extended Krylov Subspaces without Explicit Inversion
It will be shown that extended Krylov subspaces –under some assumptions– can be retrieved without any explicit inversion or system solves involved. Instead we do the necessary computations of A−1v in an implicit way using the information from an enlarged standard Krylov subspace. It is well-known that both for classical and extended Krylov spaces, direct unitary similarity transformations exist...
متن کاملShort recurrences for computing extended Krylov bases for Hermitian and unitary matrices
It is well known that the projection of a matrix A onto a Krylov subspace span { h, Ah, Ah, . . . , Ak−1h } results in a matrix of Hessenberg form. We show that the projection of the same matrix A onto an extended Krylov subspace, which is of the form span { A−krh, . . . , A−2h, A−1h,h, Ah, Ah . . . , A`h } , is a matrix of so-called extended Hessenberg form which can be characterized uniquely ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- SIAM Review
دوره 50 شماره
صفحات -
تاریخ انتشار 2008